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Abstract 

We adopt the '0-1' test for chaos using Brownian motion chains to identify the 
dynamics of the Ricker's population model. In the '0-1' test '0' is related to regular 
motion while '1' is associated with chaotic motion. The identified regular and chaotic 
types of solutions have been confirmed by means of recurrence plots. 
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The Ricker's model [1] has been invented to analyze the population of fish. 
Originally the model used for stock assessment to learn and predict the popu- 
lation of fish and to estimate how many fish can be caught without impacting 
future production of fish. The dynamics of population change is given by the 
following map [1]: 

x n+1 = F(x n ), F(x n ) = Ax n exp(-x n ), (1) 



where x n means the population of consecutive fish generations and the con- 
stant \x depends on biological conditions, F( Xfij is ci discrete transform func- 
tion. A similar model has been also used in the analysis of insect data by Moran 
[2]. This simple model give complex answer for different system parameters. 
It has attracted researchers interested in emerging the chaotic solutions in 
mathematical biology [3]. Recently, it has been also analyzed in presence of 
periodic or noisy forcing [4,5]. 

In this paper we will apply Eq. (1) to generate the time series of fish gen- 
erations which can be latter analyzed by nonlinear stochastic methods. The 
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corresponding transform functions F(x n ) for chosen values of the control pa- 
rameter fx (/i = 8, 12 18.5 and 20) are plotted in Fig. 1. 

We started examination of the dynamics of the above model with the time 
histories simulated for above chosen (Fig. 1) fj, values. The time histories are 
presented in Fig. 2a-d. First three figures Figs. 2a-c correspond to periodic 
solutions with two (Figs. 2a-b) or four points (Fig. 2c). The last figure (Fig. 2d) 
represents a chaotic solution. Note that the steady state for periodic solutions 
is reached after few initial iteration points. 
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Fig. 1. Transform function F{x n ) (Eq. 1). T' denotes the curve for [i = 8 while 
! 2'-'4'; /I = 12, 18.5 and 20, respectively. 

Naturally the standard way to classify the dynamics of nonlinear time series 
is to perform reconstruction of phase space [6] and calculate the Lyapunov 
exponent [7,8]. Alternatively one can adopt of '0-1' test method introduced by 
Gottwald and Melbourne [9]. The method was previously successfully tested 
for the Van der Pol system, the Kortweg-de Vries equations [9] the logistic 
map, and the Lorenz 96 system, including the system with weak noise [10]. 

Let us rearrange the x n series in the following way: 

n 

P( n ) = H x i cos 0' c )> ( 2 ) 

3=1 

for consecutive n = 1, 2, 3, ... . Here the constant c = 3.7 has been chosen 
arbitrary as in Ref. [9,10]. 

Note that we build new series of p(n) (Eq. 4) starting with an essentially 
bounded coordinate x n . However, the new series can be either bounded or un- 
bounded depending on dynamics of the examined process: regular or chaotic, 
respectively. In fact for a random noise T n we have a Brownian motion analogy 




2 




50 100 150 200 / H x 50 100 150 200 



n w n 

Fig. 2. Time histories of x n (a-d) for \i = 8, 12, 18.5 and 20, respectively. The initial 
point was set to x\ = 0.1. 

[11,12,13]: 

p(n + 1) = p{n) + T n n = l,2,...,N (3) 

Here p(n) represents unbounded time series where total displacement is scaled 
by the square root of steps number N [11]. T n is a diffusion term which couples 
different states n ^ m. In our case this term is given by the formula [9,10]: 

T n = x n cos(nc). (4) 

Needless to say, in case of chaotic time series x n , T n will mimic random be- 
haviour [12,13] leading to an unbounded displacement. 

On the other hand for a periodic solution x n : 

^ ^lj ^V; ^1) •••■) <^Vj ^lj ; for 5 <^ 7", (5) 
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Fig. 3. 'Phase portrait' in p — q coordinates. Here p(n),q(n) have been calculated 
for n = 1, 4000000. Arrows show the total displacement from the first to last 
points. The first point was set by assuming x\ =0.1. Note the scale in each figure 
is different. 

where the index r is related to a r-element period included in x n series. If it 
appears j times in the total time series x n then n = rj + s. The parameter 
p(n) can be expressed as (see Appendix A): 

n r+e(s-k) 

v{n) = p(jr + s) = x i cos U c ) ~ ~ 2 x k R nr sin(kc + (f){n r )), (6) 

j=l k=l 



where 4>(n r ) and R nr is the angle a characteristic circle radius (Eq. ??, here 
R nr ~ 1) dependent on the number of repetitions of a period Xx,...,x r and a 
value c (Appendix A) and 0(-) is the Heaviside step function. 

Defining the other complementary quantity q(n) in the similar way: 

n 

q( n ) = x 3 sin O) ( 7 ) 

.7 = 1 
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Fig. 4. K versus [x. Vertical straight lines correspond to fj, = 8, 12, 18.5 and 20 
respectively. 

we can derive a corresponding simplified expression for a periodic series x n 
(Eq. 5) in an analogous way. Thus an appropriate expression (Appendix A) 
reads: 



n r+e(s-i) 

q(n) = q(jr + s) = ^ Xj cos(jc) ps ^ a^i^. cos(zc + 0). (8) 
j=i i=i 



Note the pair (q,p) defines new coordinates in the two dimensional space. 
Knowing that r and s (Eqs. 6, 8) are usually low natural numbers (see Fig. 2 
r = 2 or 4) we can conclude that one can easily see if series p(n) and q(n) are 
bounded for arbitrary large n. Furthermore the simple expressions Eqs. (6) 
and (8) determine their trajectory in a circle like shape which radius R can 
be easily estimated as 

R r < Xi > . (9) 



In the above expression < ■ > denotes the average value while r is a period of 
{x n } (Eq. 5). In Fig. 3 we show examples of trajectories (plotted by consecutive 
points) for the examined system with four values of A as in Figs. 1-2 calcu- 
lated from x n series using Eqs. (2,7) and for n = 1, ...,4000000. Arrows show 
the total displacement from the first (n = 1) to last points (n = 4000000). 
Thus depending on motion type the transformed system, into p and/or q co- 
ordinates, their time series show the bounded or unbounded behaviour. 

To classify the dynamics in a more systematic way, we have performed the 
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Fig. 5. Recurrence plots of x n (a-b) for fi = 18.5 and 20 respectively. 



next step by calculation mean square displacement of p(n) [9,10]: 



N 



M{n) = lim — £ (p(j + n) - p(j)f 



(10) 



together its asymptotic behaviour of M(n) estimated by the following limit 



[9,10]: 



K = lim 



ln(M(n) + 1) 



(11) 



Inn 



Although the limits in the above two expressions go to oo, in practice, it is 
enough to assume 



to get a clear distinction between regular and chaotic states. The results for K 
versus fi are plotted in Fig. 4. It is apparent that in case of chaotic oscillations 
K 1 (mainly above /j = 14.6) while for regular motion K ps (see region 
below \i = 14.6). Some variations around and 1 can be explained by a 
finite truncation (Eq. 12). As noted by Gottwald and Melbourne [10] the 
choice of the parameter c can be sometimes important creating resonances 
which should be avoided. However we have checked this possibility in our 
calculations and assumed it eliminates such unwanted cases. Note that there 
exists a sudden change of K m 1 into lower value K m 0.38 for fi about 18.5. 
Looking at Fig. 2c we know that the system is periodic with a characteristic 
period r = 4. It is also bounded in terms of coordinates p(n) and q(n) (as 
can be seen in Fig. 3c) but the characteristic circle radius R in (p, q) space 
is relatively large (Fig. 3c, Eq. 9) comparing to other regular states (Fig. 3a- 
b). It is clear from Eq. 9 as r is relatively large here. To show the difference 
between both cases /j = 18.5 and 20 we plotted corresponding recurrence plots 
defined by the scalars p(n) (without embedding) [14]. Naturally the usual way 
is to investigate a recurrence plot in a reconstructed phase space [15]. But 
as has been recently shown by Thiel at al. [14] treatment without embedding 
gives correct estimation of dynamical invariants and various recurrence plot 
analysis measures can be applied [16]. Here we plot a dot an a squared matrix 
whose axes correspond to the % and j indices if only the condition 



is full filed, where e is a certain threshold value. Using this method we can again 
classify independently the states of the system as regular for /i = 18.5 (Fig. 
5a) with regular patterns (continuous diagonal lines) and chaotic for /i = 20, 
where the lines patterns appear with different finite length distributed in some 
random way. It indicates that the system is visiting the same region of the 
attractor many times [15]. The length of these short diagonal lines can be 
related (inversely proportional) to the largest positive Lyapunov exponent Ai 



In summary, the 0-1 test for chaos appeared to be reliable method in case of 



N = 150000 and n = 3000000 



(12) 




(13) 



[17,15]. 
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Ricker's model. Using a certain projection of the original time series {x n } into 
p(n) we could use the stochastic methods to analyze the asymptotic behaviour 
of mean square displacement (in p(n)) . It enabled us to classify the system 
dynamics directly from scalar signal without any procedure of phase space 
reconstruction. Like in other systems analyzed by Gottwald and Melbourne 
[9,10] for chaotic changes of consecutive fish populations we observed clear 
oscillations around K — 1. On the other hand in case of regular systems it is 
more useful if the period of repeating chain in a time series is shorter. However 
in longer periods (r = 4 for p, — 18.5 Fig. 4) we observed K ps 0.38. We hope 
that this disadvantage can be reduced in future by using the estimated radius 
R (Eq. 9) in the definition of K. The above method seams to be very useful 
for analysis of experimental data as well as for simulated dynamical systems 
especially those having non-continuous or non-smooth vector fields [18,19,20]. 

Acknowledgements: GL would like to thank Max Planck Institute for the 
Physics of Complex Systems in Dresden for hospitality. 

Appendix A 

For a time series represented as a large multiple sequence xi, x r 

X\, . . . , x r , X\ , . . . , x r , X\ , . . . , x s (A. 1) 

n s / n r \ r / n r \ 

Y Xj COs(jc) = Y X i[ H cos (^ c + kc) J + Y X i[ H cos (^ c + Ml ) ( A - 2 ) 
j=\ i=l \fc=l / i=s \k=l / 



where n r corresponds to the multiplicity of a shorter r-element (xi,...,x r ) 
period inside x n time series. Using a simple decoupling 

cos(ic + kc) = cos(ic) cos(/cc) — sin(ic) sm(kc) (A. 3) 



and 



A nr = cos(kc) (A.4) 

k=l 

cose — cosc(?7, r + 1) — 1 + cos c cos c(n r + 1) + sin c sin c{n T + 1) 

2 — 2 cos c 



Similarly 



B nr = Yl sin (M 



(A.5) 



k=l 



— sin c cos c(n r + 1) + cos c sin c{n r + 1) + sin c 
2 — 2 cos c 
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Figure A.l presents results B nr versus A nr for two different values of c (c = 1.7 
and c = 3.7) note that in both cases the points are lying on circles. They are 
however centered in different points and have different radii R c . In case of 
c = 3.7 the radius is less than f . The centers are given by the averages taken 
from Eq. A.4 and A. 5. 

A nr o = 0.5 for any c, (A. 6) 

while 

sin c 

B nr = ~Z ~ • (A. 7) 

2 — 2 cos c 

Note, using the circles representation, both expressions for A nr and B nr (Eqs. 
A.4 and A. 5, respectively) can be now written in a compact form 

A nr = A nr0 + R c cosip(n r ), 

B nr = B nr0 + R c smif)(n r ). (A. 8) 

Interestingly for c = 3.7 we have found numerically a change of if) versus n r is 
linear 

ip(n r ) « 0.5584n r . (A.9) 
Using the above notation we can write 



„ r+G(s-k) 

Xj cos(jc) = x k (A nr cos(ck) - B rn sin(fcc)) 

j=l k=l 
r+Q(s-k) 

= ^2 x kRnr COs(c/c + 4>{n r )) (A. 10) 

k=l 



where 



tan(0) = ^, R nr = jAj r + B^ (A.11) 

JD nr 



where G(-) is a Heaviside step function. 
Similar analysis can be applied for 



r / n r 



Y Xj sin(jc) = Y x i\ Y sin (^ c + kc) \ + x i [Yl sin (^ c + kc) j , (A. 12) 

j=l i=l \k=l ) i=s \k=l ) 
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Fig. A.l. B nr versus A nr for c = 1.7 and c = 3.7, respectively. Arrows indicate 
corresponding circle centers. 

where 

sin(ic + kc) = sin(zc) cos(A;c) + cos(zc) sin(fcc). (A. 13) 

Therefore 

n r+e(s-k) 

Xj cos(jc) = x k (A nr sm(ck) + B rn sin (fee)) 

3=1 k=l 
r+&(s-k) 

= Y, x kRnr sin(c/c + (f)(n r )). (A. 14) 

k=l 
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